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Abstract. Coude spectroscopic data of a poorly-studied 
peculiar supergiant, HD327083, are presented. Ha and H/3 
line profiles have been fitted employing a non-LTE code 
adequate for spherically expanding atmospheres. Line fits 
lead to estimates of physical parameters. These parame- 
ters suggest that HD327083 may be close to the Luminous 
Blue Variable phase but it is also possible that it could be 
a B[e] Supergiant. 
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1. Introduction 

Luminous emission-line B-type stars often present mas- 
sive winds. Many of them share some common proper- 
ties: i) optical spectra dominated by HI and Fell emis- 
sion lines, usually displaying P-Cygni profiles; ii) presence 
of forbidden emission lines, specially of [Fell]; iii) strong 
red/infrared excess. Such properties are typical of the B[e] 
Supergiants (Lamers et al. 1998) but are also often seen in 
the Luminous Blue Variables (LBVs). Their localization 
in the H-R diagram is a real problem. The distances are 
very uncertain (particularly for galactic objects) and the 
absence of photosphcric lines renders the assignment of a 
spectral type quite difficult. 

The present paper concerns one of such interesting 
stars, namely HD327083. Almost no study exists in the lit- 
erature devoted to it, but some works on emission- line ob- 
jects have included it. Ha emission was discovered by Mer- 
rill & Burwell (1949). Carlson & Henize (1979) included 
HD327083 in their list of peculiar southern emission-line 
objects. They have proposed a B8 spectral class. Some 
years ago Lopes et al. (1992) showed observations of some 
spectral regions. They revealed HI lines with P-Cygni pro- 
files as well as Fell AA 4924, 5018, 5169 A (m42). On the 
other hand the Fell A transition at 9997 A seemed to be 
in pure emission. The authors have suggested a B6 type. 
More recently. The et al. (1994) presented a new cata- 
logue of (confirmed and candidates) Hcrbig AeBe stars. 
In this catalogue HD327083 is cited among the extreme 
emission-line objects, which are not identified to belong 
to a certain group. However, they proposed the spectral 
type B1.5. The same assignment has been used by Valenti 
et al. (2000) while Sheikina et al. (1999) have obtained log 
Teff > 4.2 and log - 6.0. 

Here we show our recent spectroscopic data and es- 
timates of some physical parameters (as M, i*, and 
^He, the numerical abimdance ratio n(He)/n(H)). Ha and 
H/3 line profiles have been fitted with this aim. These es- 
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Fig. 1. Normahzed spectra: Ac = 4750, 5000 A (see text). 



timates allow us to locate (within certain uncertainties) 
HD327083 in the H-R diagram. A discussion of its evo- 
lutionary status, based on the derived parameters, is also 
presented. 



2. The data 

The spectra were obtained in August, 1999 using a Coude 
spectrograph at the LNA 1.6m telescope in Pico dos 
Dias (Brazil). A grating of 600 1/mm combined with a 
1024 X 1024 CCD have been used; the reciprocal dispersion 
is about 0.4 A pixel"^ yielding a wavelength range of ~ 
400 A . The central wavelengths are: Ac = 4750, 5000, 5200 
and 6400 A. The data have been reduced using IRAF rou- 
tines. We followed standard procedures: bias subtraction, 
flat-fielding, linearization, wavelength calibration and cor- 
rection for atmospheric extinction. 
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Table 1. Fell lines 
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Fig. 2. Normalized spectra: Ac = 5200, 6400 A (see text). 

In Figs. 1-2 we can see our normalized spectra. In fact 
the spectral regions covered by our observations are domi- 
nated - besides Ha and H/3 -, by Fell emission lines. Some 
of these lines show a P-Cygni profile and others are in 
pure emission. Table 1 lists the Fell transitions identified, 
which arc marked in Figs. 1-2. Column 1 gives the mea- 
sured wavelength (A) of emission peaks, followed by our 
identification (Colmim 2). Wc note that only one feature, 
at A 5159.0 A, has been recognized as a forbidden transi- 
tion (see Fig. lb). On the other hand [01] A 6300.2 A is 
present (Fig. 2b). 

3. Summary of the model 

Let us briefly describe now the model employed. Re- 
cently Machado (1998, 1999) developed a non-LTE nu- 
merical code for radiation transfer in expanding stellar 
atmospheres of massive stars. This code is based on basic 
assumptions: spherical symmetry, stationarity and homo- 
geneity. The density structure is related to the mass loss 
rate and the velocity field V{r) via the equation of conti- 
nuity 

= 

with M as an input parameter and the velocity field is 
pre-specified in an ad-hoc way as a /3-type law 
V{r) = V^{l-RWrf (2) 
where Voo is the terminal velocity of the wind and /3 is a 
measure of the rate of acceleration (de Koter et al. 1993). 

The code is based on the following scheme. The sta- 
tistical equilibrium equations are solved using the escape 
probability method for calculating the source function 
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while the transfer equation is solved using the "SEI" - 
Sobolev Exact Integration - method (Lamers et al.l987). 
This is justified by the flow velocities observed in the wind 
of HD327083. 

In the present version it is considered that the wind 
begins in the sonic point localized at radius i?*. The effec- 
tive temperature is defined by the luminosity and 
the radius R^, via the Stefan-Boltzmann law. 

Regarding the temperature structure we have used the 
law suggested by Drew (1985): 

v{ry 



T{r) 



0.78- 



0.51- 



(3) 



We are aware that this temperature profile is more appro- 
priate to main sequence O stars. Thus we have checked 
the influence of the temperature law. Some simple tests 
with isothermal and radially decreasing ad-hoc profiles 
have been performed. It was found that the final result 
is not sensitive to similar temperature laws. 

The model atmosphere in the present work consists of 
hydrogen and helium. For H and Hell we consider the low- 
est ten bound levels. For Hel we adopt the atomic model 
suggested by Almog & Netzer (1989) that includes the in- 
dividual levels until n <4 and combined levels with 5 < n 
< 10. 
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The determination of the population numbers of the 
energy levels is made in a non-LTE way solving the equa- 
tions of statistical equilibrium for each level considered 
in the atomic model. The radiative transition probabilities 
and the photoionization cross section for all levels consid- 
ered are quoted from the Opacity Project. 

The ionization collision cross sections have been ob- 
tained from the approximate formulae given by Jefferies 
(1968). Excitation collisional rates have been calculated 
from the collisional strength values. For H and Hell these 
values are computed using polynomials given by Giova- 
nardi et al. (1987). For Hel the collision strengths are 
taken from Bcrrington & Kingston (1987) for levels with 
n <4 and from Auer & Mihalas (1968) and Mihalas & 
Stone (1968) to the other ones. 

The radiation and the statistical equation systems are 
solved simultaneously through the wind that is divided in 
shells. A convergence criterium is applied in order to move 
from one shell to the next. 
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4. Ha and Hj3 line fits 

4-.1. Methodology 

The input parameters are the temperature T^, the stellar 
radius i?* (and luminosity), the mass loss rate M, the 
terminal velocity Voo , the helium abundance Ahc and the 
free parameter of the velocity profile /3. In the present 
work we have fixed the value of the terminal velocity, l^o 
= 400 kms~^. This value is an average of the Vedges of Ha 
and H/3 profiles. 

Due to the great number of "free" parameters we have 
adopted the following strategy in order to obtain the line 
fits. First we defined temperature limits. Previous deter- 
minations of temperature for HD327083 are very uncer- 
tain. As discussed in the introduction first estimates have 
pointed to values around T^b ~ llOOOK - 13000K. More 
recent ones have lead to larger values: Tofi ^ 16000K - 
19000K. So we have decided to explore a broad range, that 
roughly corresponds to the LBVs and Yellow Supergiants 
domain: 8000K - 31000K. This Toff range was fixed based 
on the works by Gummersbach et al. (1995), Szeifert et 
al. (1996) and de Jager & Nieuwenhuijzcn (1997). 

Secondly we have started computing models according 
to the points of the evolutionary tracks by Schaller et al. 
(1992) from Mzams = 20 M© to 120 M©. Each track con- 
sists of a number of points which provide stellar temper- 
ature, luminosity, mass loss rate and surface abundance 
ratio (^Hc)- These arc the entry parameters for the code. 
The evolutionary points which give reasonably good fits 
(according to a criterium) are picked as starting points. 
From them we subsequently varied the input parameters, 
specially mass loss, in order to obtain better adjustments. 



Fig. 3. Observational data (dashed lines) and best theo- 
retical fits (full lines) for Ha and H/3 profiles. 

4-2. Results 

The mass loss rate is indeed small for the Mzams = 20 
Mq and Mzams = 25 Mq paths. It is not enough to pro- 
duce Ha and H/3 profiles with the intensities observed 
in HD327083. Concerning the Mzams = 85 M© and the 
Mzams = 120 M© tracks, in such cases the He abundance 
is already high for stars just evolved off the Main Se- 
quence. In these regions the temperatures are high (Ti, > 
40000K) and appropriate fits of Ha and H/3 cannot be 
achieved with reasonable values for the others parameters. 

Our best fits have been obtained with the following set 
of parameters. Model 1: M = 4.9 x 10"^ M© yr'^ = 
19000K, L^ = 5x 10^ L©, ^He = 0.40. Model 2: M = 1.0 
X 10-'' M© yr-i, = 9000K, = 9.5 x 10^ i©, Ahc = 
0.11. These fits are shown in Figs. 3a and 3b respectively 
and were both obtained with /3 = 3, which indicates a 
slow expansion. Model 1 is based on the series of Ha and 
H/3 profiles along the Mzams = 40 M© path while model 
2 is based on the Mzams = 60 M© path. They have been 
selected following the methodology described earlier. We 
note that the "best model" values of M are slightly higher 
(about 50%) than those given directly by the evolutionary 
paths. 

It can be seen that the agreement between observa- 
tional and theoretical profiles is rather poor. The models 
reproduce well the emission strengths but the absorption 
features are quite flat and narrow. This could suggest that 
the emission originates in a high density disk. Unfortu- 
nately we cannot check this alternative within the scope 
of our spherical model. 
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Fig. 4. Theoretical evolutionary tracks of 25 M©; 40 M©; 
60 Mr;: and 85 M0 initial main-sequence mass. The circles 
locate the best model parameters. 



5. Discussion 

Each set of parameters leads to a location in the H-R dia- 
gram. Fig. 4 shows these positions and some evolutionary 
paths. Two alternatives - redward or blueward - are pos- 
sible if only and are considered. He abundance may 
be invoked to choose one possibility or another. 

The evolutionary model 1 has Ahc = 0.40. This value 
is in the blueward track: the star would be coming back to 
the hotter region of the diagram. The parameters in this 
scenario point to a LBV quiescent (or post-LBV) phase. 
On the other hand the evolutionary model 2 has j4hc = 
0.11. This low value indicates a point in the redward track. 
Here the star is in a less evolved stage, coming to the LBV 
eruptive phase. The intensity of the He lines should, in 
principle, reflect the He abundance. He lines are not seen 
in our data. We note however that our spectra do not 
cover the wavelength regions with strong He lines, such as 
A = 5876, 6678 A. So, we prefer not to discard model 1 at 
present time. 

Both alternatives allow us to propose that HD327083 

is near the LBV phase. However, this suggestion must be 
taken with caution. First because we have fitted only Ha 
and H/3 line profiles. So, the uncertainty in the derived 
parameters is large. Second because the high value of (3 
needed to obtain "good" fits reinforces the scenario of 
emission coming mostly from a flattened region, as the 
equatorial disk believed to be present in B[e] Supergiants 
(Zickgraf et al. 1986). 



In general the B[e] stars have stronger [FeH] lines than 
the LBVs. In our spectra only the feature at A 5159.0 A 
(which we attribute to transitions in multiplets [18F] and 
[19F]) seems to be present in Fig. lb, although it is im- 
certain in Fig. 2a. It is true that our data only covers the 
regions A ~ 4550 - 5400 A and A ~ 6200 - 6000 A. How- 
ever, several forbidden transitions (from multiplets [3F], 
[4F], [18F], [19F], [35F] etc) are in these regions. On the 
other hand [01] A 6300.2 A is clearly seen (Fig. 2b). This 
feature is intense in the B[e] Supergiants but absent in the 
LBVs. 

It is likely that HD327083 is in a short-lived evolved 

stage of massive stars. However it is still unclear if it is a 
LBV candidate or if it is a B[e] Supergiant. It is sure that 
it deserves to be further investigated. Such analysis should 
be based on higher-quality spectroscopic data, covering a 
broad spectral domain. Moreover it should include other 
information about the star such as the energy distribution, 
Balmer jump, and colour. From the theoretical point of 
view it is needed to employ a non-spherical code(e.g. Stee 
et al. 1995) to obtain adequate fittings not only for HI lines 
but also for Hel and N lines. This detailed investigation 
will clarify the nature of HD327083 and it may contribute 
to a better understanding of evolved massive stars. 
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